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ABSTRACT 

Aims. We explore the accretion disk and boundary layer emission in bright neutron star LMXBs and their dependence on the 
mass accretion rate. 

Methods. Fourier-frequency resolved spectroscopy of the archival RXTE data. 

Results. We demonstrate that Fourier-frequency resolved spectra of atoll and Z- sources are identical, despite significant 
difference in their average spectra and luminosity (by a factor of ~ 10 — 20). This result fits in the picture we suggested earlier, 
namely that the / ^ 1 Hz variability in luminous LMXBs is primarily due to variations of the boundary layer luminosity. In 
this picture the frequency resolved spectrum equals the boundary layer spectrum, which therefore can be straightforwardly 
determnined from the data. The obtained so boundary layer spectrum is well approximated by the saturated Comptonization 
model, its high energy cut-off follows kT ~ 2.4 keV black body. Its independence on the global mass accretion rate lends 
support to the theoretical suggestion by Inogamov&Sunyaev (1999) that the boundary layer is radiation pressure supported. 
With this assumption we constrain the gravity on the neutron star surface and its mass and radius. Equipped with the 
knowledge of the boundary layer spectrum we attempt to relate the motion along the Z-track to changes of physically 
meaningful parameters. Our results suggest that the contribution of the boundary layer to the observed emission decreases 
along the Z-track from conventional ~ 50% on the horizontal branch to a rather small number on the normal branch. 
This decrease can be caused, for example, by obscuration of the boundary layer by the geometrically thick accretion disk 
at M ~ A^Edd- Alternatively, this can indicate significant change of the structure of the accretion fiow at M ~ Msdd and 
disappearance of the boundary layer as a distinct region of the significant energy release associated with the neutron star surface. 
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1. Introduction 

Accreting neutron stars in low mass X-ray binaries 
(LMXB) are among the most luminous compact X-ray 
sources in the Galaxy. At least several of them have lu- 
minosities exceeding ~ few x 10'^* erg/s and presumably 
accrete matter at the level close to the critical Eddington 
accretion rate. E arly observations of these sources (e.g. 
iToor et aliFlOTOl) revealed rather soft X-ray spectra, indi- 
cating that their X-ray emission is predominantly formed 
in the optically thick media. Similar to accreting black 
holes, at lower X-ray luminosities (lower mass accretion 
rates), Lx ^ 5 x lO^^^ erg/s, neutron st ars undergo a tran- 
sition to the hard spectral state (e.g. lBarretll200l|l . The 
energy spectra in this state point at the low optical depth 
in the emission region. 

In the soft spectral state, the commonly accepted 
picture of accretion at not too extreme values of ac- 
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cretion rate has two main ingredients - the accre- 
tion disk (AD) and the boundary layer (BL). While 
in the disk matter rotates with nearly Keplerian ve- 
locities, in the boundary layer it decelerates down to 
the spin frequency of the neutron star and settles 
onto its surface. For the typical neutron star spin fre- 
quency 500 — 700Hz) comparable amounts of energy 
are released in these two region s (^Sunyacy & Shakura^ 
Il986t ISibgatullin fc Sunvaevll2000|) . This picture is based 
on rather obvious qualitative expectations as well as 
more sophistic ated theoretical considerat i ons and numer- 



ical modeling JSunvaev & Shakuralll986l; iKluzniak 


|l988l 


llnogamov & Sunvaevlll999l ISibgatullin & SunvaevI 


2ood). 



It has been receiving, however, little direct observational 
confirmation. Due to similarity of the spectra of the ac- 
cretion disk and boundary layer the total spectrum has 
a smooth curved shape, which is d ifBcult to decoinpose 
into separate spectral componeiits JMitsuda et alj [19841: 
IWhite et alJll98aHDi Salvo et al.ll200lUDone et alJl2002l) . 
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This made application of physically motivated spectral 
models to the description of observed spectra of luminous 
neutron stars difficult, in spite of very significant increase 
in the sensitivity of X-ray instrume nts. A possibl e solution 
was suggested by early results of iMitsuda et al.l l|l984f) . 
which demonstrated the potential of using the combined 

spectral and variability information. 

Recently, ICilfanov. Revnivtsev & Molkovl l(200.ll) ana- 
lyzed spectral variability in luminous LMXBs and showed 
that in these sources aperiodic and quasi periodic vari- 
ability on ~ sec - msec time scales is caused primarily 
by variations of the luminosity of the boundary layer. Its 
spectral shape remains nearly constant in the course of 
the luminosity variations and is represented by the Fourier 
frequency resolved spectrum. Moreover, in the considered 
range M '-^ (0.1 — l)MEdd {MEdd is the critical Eddington 
mass accretion rate) it depends weakly on the global mass 
accretion rate and in the limit M ^ Mscid is close to Wien 
spectrum with kT ~ 2.4 keV. Such a behavior is in accord 
with the predictions of the model by Inogamov & Sunyaev 
(1999), namely, that at sufficiently high accretion rates, 
M > O.lMEdd, the boundary layer is radiation pressure 
dominated and the local radiation ffux is close to the crit- 
ical Eddington value. Increase of the mass accretion rate 
leads to the increase of the emittin g area of the BL, while 
its ve rtical structure changes little ijinogamov fc Sunvaevl 
Il999t) . 

In this paper we further explore the behavior of the 
boundary layer and the accretion disk emission in a num- 
ber of bright neutron star LMXBs (both Z and atoll 
sources) and qualitatively describe their evolution with 
change of a mass accretion rate. Our study will b e base d 
on the results of ICilfanov. R.evnivtsev fc Molkovl ((200,ll) . 
namely, that the shape of the boundary layer spectrum 
is adequately represented by the frequency resolved spec- 
trum (energy spectrum of the variable part of the emis- 
sion) at the Fourier frequencies above / ^ 1 Hz. We 
will concentrate on neutron star LMXBs in the soft/high 
spectral state with luminosities Lx 0.5 — 1 x 10'^'' 
erg/s. Sources in this state have very impotant property - 
their frequency resolved energy spectra do not depend on 
Fourier frequency. It is this property which allow us to de- 
compose the averaged energy spectrum of the sources into 
two components. In hard spectral state the variability is 
significantl y more comphcated (see e.g. lR,evnivtsev et all 
Il999il200(ll) and we leave it for other studies. 

2. Data 

For our study w e used data of RXTE observatory 
ifBradt et a,1.lll 99r!|> . which combines large collecting area 
with high time resolution. We have studied all Z-sources 
except for GX349-I-2 which during the majority of RXTE 
observations was on the "ffaring" branch of its color-color 
diagram, supposedly corresponding to super- Eddi ngton 
mass accretion rates llHasinger fc van der KliFj Fmsl. The 
behavior of sources at such high accretion rates is more 
complex and is beyond the scope of this paper. The re- 



maining five Z sources are: Cyg X-2, GX340H-0, Sco X-1, 
GX 5-1 and GX 17-f 2. For these sources we have excluded 
the time intervals corresponding to the "fiaring" branch 
of their color-color diagram. The atoll sources were repre- 
sented in our sample by 4U1608-52 and 4U1820-30 in their 
soft/high spectral state, showing high rapid flux variabil- 
ity. 

For the data reduction we used standard programs 
of package HEASOFT 5.3 in accordance to the RXTE 
Guest Observer Facility recommendations. We have cor- 
rected the obtained energy spectra for the deadtime 
and the pile up effects ( "http:/ /lheawww.gsfc.nasa.gov/( 
docs/xray/xte/pca/). All subsequent spectral approx- 
imations included interstellar absorption which was 

6 X 10^2 cm-2 for GX 

,-2 



fixed a t values : Nh 

M O+0 lllaria et all l2004|) . 0.2 x 10^^ cm^^ for Cyg X- 
2 teuulkers et alJ Il997l). 0.2 x lO^^ cm^^ for Sco X- 
1 l|Kahn. Sewa rd, fc Chlebowski"l984V 1.0 x 10^^ cm-^ 
for GX 17+2 (Langmcicr. Hasingcr. fc TruemDerlll99nl) . 
3.0 X 10^2 cm-2 for G X 5-1 (Asai et allll994ll . 1.0 x 10^^ 
cm-2 for 4U1608-52 (Penninx et al."l989') and 0.2 x 10^^ 
cm-2 for 4U1820-30 ( Habcrl et a,l. 1987). The luminosi- 
ties of all Z sources are approximately 2 x 10'^^ erg s~^, 
for 4U1608-52 and 4U1820-30 approximately lO^"^ erg s'^ 
The color-color diagrams (CCD) of these sources are 
presented in Fig^ For construction of the CCDs we have 
used energy channels 2.0-4.0 keV, 4.0-6.0 keV, 6.0-10.0 
keV and 10.0-15.0 keV. The hard and soft colors are ra- 
tio of energy fluxes of sources in abovementioned energy 
bands. Energy fiuxes were computed from spectra aver- 
aged over 128 sec intervals. All fluxes were corrected for 
absorption. 



3. Spectrum of the boundary layer 

The frequency resolved spectra ( the energy spectrum o f 
the variable part of the emission, iR.evnivtsev et al]ll999l) 
computed at the frequencies / ^ 1 — 5 Hz are shown 
for five sources from our sample in Fig|21 For Z sources 
we used only data on the horizontal branch of the CCD 
where the amplitude of variability at these frequencies is 
maximal. For Sco X-1 and GX 5-1 the RXTE achive does 
not contain high time resolution data for the horizontal 
branch with sufficient number of energy channels. The fre- 
quency resolved spectrum of 4U1608- 52 was taken from 
ICilfanov. Revnivtsev fc Molkovl l)2003|) . For plotting pur- 
poses, the normalizations of all spectra were adjusted to 
match that of GX340+0. 

Similarity of the spectra shown in Fig|21is remarkable, 
considering significant difference in the average spectra 
and a factor of ^ 10 — 20 spread in the luminosity be- 
tween atoll and Z-sources (~ O.lMsdd and ~ Msdd corre- 
spondingly). This behavior fits in the p icture proposed by 
ICilfanov. R.evnivtsev fc Molkovl l|2003l) . As summarized in 
the Introduction, they showed, in particular, that the fre- 
quency resolved spectra at these frequencies equal the 
spectrum of the boundary layer emission. 
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Fig. 1. Color-color diagram of all Z-sources except for 
GX349-I-2. Only horizontal and normal branch data are 
shown. Colors were calculated as ratios of absorption cor- 
rected energy fluxes in the following energy channels 3.0- 
4.0 keV, 4.0-6.0 keV, 6.0-10.0 keV, 10-15.0 keV. Dashed 
line shows the approximation to the Z-track by an ar- 
bitrary smooth curve. Turning points of the Z-track are 
marked with Sz = 0, 1, 2. Here and later black points de- 
note GX340-hO, red - Cyg X-2, cyan - GX 5-1, magenta - 
Sco X-1, blue - GX17-I-2 

The shape of the frequency resolved boundary 
layer) spectrum can be adequately described by the sat- 
urated Comptonization. For the sake of comparison with 
other results and for convenient parametrization of the BL 
spectrum we used the Comptonization model of Titarchuk 
(1994). The best fit parameters of the model fitted in the 
3-20 keV range simultaneously to all five spectra shown in 
the FigElare: temperature of seed photons kTg ~ 1.5±0.1, 
temperature of electrons kT^ = 3.3 ± 0.4 and the optical 
depth r = 5 ± 1 for slab geometry. The best fit model is 
shown by the thick dotted line on Fig[21 The temperature 
of the black body spectrum describing the shape of the 
cutoff in the observed spectrum at energies >13 keV is 
fcTbb = 2.4 ± 0.1 keV (thin dashed line on Figl^Jl. 

It would be interesting to fol l ow up on the results 
of iGilfanov. R,evnivtsev & Molkovl ( 20031) and to consider 
the evolution of the boundary spectrum along the Z track 
from the horizontal to normal branch. Unfortunately, 
apart from GX340+0, already considered in their paper, 
the variability level on the normal branch of color-color 
diagram for other four Z-sources is insufficient to obtain 
frequency resolved spectra with reasonable signal-to-noise 
ratio. 



Fig. 2. Fourier-frequency resolved spectra (^boundary 
layer spectra) of sources from our sample. All spectra were 
corrected for the interstellar absorption. The thick dot- 
ted line shows the best fit Comptonization model with 
kTs = 1.5, kTe — 3.3 keV, r = 5. Thin dashed line shows 
the blackbody spectrum with temperature /cTbb — 2.4 
keV. 



4. Constrains on NS mass and radius from the BL 
spectrum 

The independence of the spectrum of the boundary layer 
on the l uminosity lends suppor t to th e theoretical predic- 
tions by llnoeamov fc SunvaevI l)l999^ that the boundary 
layer is radiation pressure supported. It means that in the 
upper layers of the BL only the radiation pressure coun- 
teracts the gravitational force. Therefore we can say that 
every unit area of the BL should emit Eddington surface 
flux. At smaller mass accretion rates (M ~ O.lMEdd) the 
BL occupies only a part of the NS surface, giving rise to 
moderate total X-ray luminosities (Lx ~ 10'^^ erg s~^). If 
mass accretion rate increases the total BL luminosity ap- 
proaches the Eddington limit for the NS iEdd ~ 2 x 10'^^ 
erg s~^ (1.4M0)^^. In this picture, the parameters of the 
BL emission can be used to determine the value of the 
Eddington flux limit on the surface of the neutron star. 
As the Eddington flux limit is uniquely determined by the 
neutron star surface gravity and the atmospheric chemi- 
cal composition, the neutron star mass and radius can be 
constrained. 

If the boundary layer emitted true black body emis- 
sion, the radiation flux of the unit area was determined 
only by its temperature. Therefore the observed shape 
of the BL spectrum, in particular the best fit black 
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by introducing the spectral hardening factor which relates 
the color and the effe c tive t e mperatures of the emission 



Neutron star radius, km 



Fig. 3. Constrains on the mass and radius of the neutron 
star in bright LMXBs. Solid lines denote the regions where 
neutron star radius R equals to 1.5i?g (radius of photon 
orbit in Schwarzschild metric) and 3.0i?g (radius of last 
stable orbit of an accretion disk in Schwarzschild metric) 



body temperature, could be used to determine the value 
of the Eddington flux limit on the neutron star sur- 
face. This approach has been utilized in th e context 
of the Eddington limited X-ray bursts (e.g. iGoldmanI 
19791 lMarshal]|ll982l: iLewin. van Paradiis. &: Taamlll993l: 



Titarchuk fc Shaposhnikovl2002() . For fully ionized hydro- 
gen atmosphere: 



aT* (TT GMil - Rsch/R) 



3/2 



i?2 



where a - Stefan-Boltzmann constant, ctt - Thomson 
cross-section, T - black body temperature at infinity, rup 
- proton mass, M - mass and R - radius of the neu- 
tron star, i?sch = 2GM/c^ - Schwarzschild radius of the 
neutron star. In addition, one would have to take into ac- 
count that the value of the Eddington flux is somewhat 
(by ^ 10 — 20%) reduced because of the action of the cen- 
trifugal force caused by the r otation of the boundary layer 
ijlnogamov fc Sunvaevll99 ^ . Rotation of the neutron star 
at this point is not very important unless it is very high 
(rotational frequency <;800-1000 Hz) 

In reality, scatterings are important in the 
atmosphere of the neutron star, therefore, the 
boundary layer sp ectrum will differ from the 
black body (e.g. iLo ndon. H oward, fc TaamI [l986; 
iLewin. van Paradiis. fc Taam 1 99,C The radiation trans- 
fer problem in the atmosphere of the neutron star has 
been intensively investigated, in particular in the context 
of X-ray bursts. Numerical calculations show that the 
effects of scatterings can be approximately accounted for 



llSvunvaev fc Shakural Il974 iLondon. H oward, fc TaamI 
1986t iTitarchukI Il994t IShimura fcTakah ara 1995; 



Ross fc FabianI I 1996(l . Typical values of the hardenmg 
factor are about ~ 1.7. We used this result in order to 
make simple estimates of the gravity on the neutron star 
surface and to constrain its mass and radius. In these 
calculations we assumed the color temperature of the 
boundary layer emission T = 2.4 keV and considered 
the range of the hardening factor values of 1.6-1.8. 
We assumed that the centrifugal force reduces the 
Eddington flux limit in comparison with non-rotating 
boundary layer by 20%. We also took into account the 
finite h eight of the boundary /spre ading layer, about 
1 km ijinogamov fc SunvaevI Il999|) . The calculations 
were performed for Schwarzschild geometry, hydrogen 
atmosphere. The result is shown in FigOl The width of 
the shaded region is defined by the assumed range of 
the values of the hardening factor. Detailed modeling of 
the emergent spectrum of the BL will allow to strongly 
diminish the size of uncertainty region. 

Finally we note that the similarity of the spectral shape 
of the BL spectrum in different sources (FiglSJ indicates 
that they have close values of the mass and radius, in par- 
ticular, that there is no significant difference in the surface 
gravity between atoll and Z-sources. It also shows that 
there are no significant differences caused by variations 
in the atmospheric chemical abundances between sources. 
In particular we did not find statistically significant dif- 
ference between ultracompact compact binary 4U1820-30 
and other sources. 



5. Evolution along Z-track 

In order to study the behavior of Z-sources in the color- 
color diagram we consider their spectra integrated over 
128-sec time intervals. These spectra are fitted with a sim- 
ple model consisting of two components, representing con- 
tributions of the boundary layer and the accretion disk. 
Also included in the model was the gaussian line at the 
energy 6.4 kcV. It is not energetically important but sig- 
nificantly improves quality of the fit. The line parame- 
ters will not be considered here. The shape of the bound- 
ary layer spectrum was approximated by comptt model 
with the best fit parameters obtained in the section 13 
For the accreti on disk spectrum we adop t the multicolor 
disk b lackbodv lshak^ira fc SunvaevI ( 197.ll) ; iMitsuda et alJ 
l|l984l) {diskhh model in XSPEC fitting package). 

Quality of the spectral fits is uniformly good on the 
horizontal and normal branches of the Z-track in the color- 
color diagram, with the reduced never exceeding ~1.3 
for 41 degrees of freedom. Relative deviations of the model 
from the data does not exceed approximately one percent. 
This is illustrated by Fig0I showing two spectra of Cyg 
X-2, at the beginning of the horizontal branch and at the 
end of the normal branch. 
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Energy, keV 

Fig. 4. Spectra of Cyg X-2 at the beginning of horizontal 
branch [Sz ~ 0) and at the end of the normal branch 
{Sz ~ 2). Solid lines are best fit models of the spectra. 
Dashed lines denote the boundary layer components 



The model is obviously oversimplified. Firstly, it re- 
lies on the assumption about the shape of the boundary 
layer spectrum. This seems to have been established rather 
well for the horizontal branch and for the upper part of 
the normal branch. For the latter, it has been shown for 
GX340-I-0 only, for other sources the signal-to-noise ra- 
tio being insufficient. On the lower part of the normal 
branch the variability level and, consequently the signal- 
to-noise ratio of the frequency resolved spectra, are too low 
to conduct a similar study for any of the sources in our 
sample. Secondly, we chose the diskbb model to represent 
the spectrum of the accretion disk. Although the model 
involves several well known simplifying assumptions it has 
been shown to approximately reproduce the shape of the 
spectrum of the geometrically thin, optically thick accre- 
tion disk (Merloni, Fabian & Ross, 2000). The best fit 
color temperature obtained in this model can be related 
with certain accuracy to the maximum disk temperature. 
Therefore with this model we should be able to separate 
the disk and boundary layer components and to estimate 
their respective fluxes, as well as the characteristic temper- 
ature in the inner disk. Keeping in mind the limitations 
of the spectral model we will attempt to relate the mo- 
tion of a source along the Z-track to the changes in the 
values of the physically meaningful parameters. The most 
important difference of our approach from similar stud- 
ies undertaken previously is in a priori knowledge of the 
shape of the boimdary layer spectrum. 




10 20 

Disk flux (0.01-30 keV) 



Fig. 5. Dependence of the disk temperature on its flux 
for five Z-sources. The units of flux are 10~^ erg/s/cm^. 
Regions corresponding to the horizontal and normal 
branches are marked. Dashed curved line shows T oc F^/J^ 
dependence. For GX5-1, GX340-hO, GX 17+2 and Sco X-1 
the flux values were scaled down by the factor of 1.3, 1.4, 
1.4 and 13 respectively to match Cyg X-2 points at the 
end of the horizontal branch 



5.1. Accretion disk 

The relation between the accretion disk temperature and 
its flux for all five Z-sources is shown in Fig|5l (we em- 
phasize that plotted along the x-axis is the disk flux, as 
opposite to the total flux). For three of them (GX 5-1, Cyg 
X-2 and GX 340+0) the disk temperature on the horizon- 
tal branch approximately follows the T oc F^-^^^ law, which 
is to be expected for a standard accretion disk. For Sco 
X-1 and GX17-f 2 the temperature increases with the disk 
flux faster than the T oc -F^isk ^'^^ sources the 
monotonic relation between the disk temperature and the 
disk flux breaks down on the normal branch and further 
along the Z-track the temperature either remains constant 
or decreases with the disk flux. 

The significant deviation from T oc F^^^ law near the 
transition from the horizontal to the normal branch indi- 
cates that the simple accretion disk model in the form rep- 
resented by xspec's diskbb becomes inapplicable. The ex- 
act reason for this is not clear. It may happen due to rather 
small modifications to the model, e.g. violation of the as- 
sumption that the effect of Comptonization can be repre- 
sented by the single value of the spectral hardening factor, 
constant over the extend of the inner disk as well as along 
the Z-track in the color-color diagram. Alternatively, this 
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Fig. 6. Dependence of the boundary layer contribution to 
the total X-ray emission of Z sources as a faction of the 
position on their Z-diagrams 

might be an indication of the more significant break down 
of the standard geometrically thin, optically thick picture 
of the accretion disk. This would point at the significant 
modification of the structure and geometry of the accre- 
tion flow near the transition from the horizontal to the 
normal branch, presumably corresponding to Ad ^ AfEdd- 
Knowledge of the boundary layer and the accretion 
disk spectra allows one to (crudely) estimate the scale- 
height of the disk H/R in the region of the main energy 
release. Indeed, for the radiation pressure supported ac- 
cretion disk its height is determined by the equality of 
the radiation pressure and the gravitational force (see e.g. 
IShakura fc SunvaevlllQySl) . For geometrically thin accre- 
tion disk the gravitational force near the neutron star 
is ~ R/H times smaller than on the neutron star sur- 
face. Assuming that the radiation transfer conditions are 
similar in the disk and in the boundary layer (i.e. the 
color hardening factors are similar, see section^ we find 
H/R ~ {Tdisk/TBh)* ~ 0.1-0.2 which is comparable with 
the predictions of the siniple dis k theory for M ^ Msdd 
((Shakura fc Sunvaevlll97ill97(il) . On the other hand, this 
estimate can be considered as a crude consistency check 
for the obtained values of the boundary layer and disk 
temperatures. 

5.2. Boundary layer 

The dependences of the BL contribution to the total X- 
ray emission on the position on the Z-track are plotted in 
FigEl The coordinate along the Z-track was defined to be 
proportional to the hard color with the reference points 
Sz — 0, 1, 2 defined as shown in FigQ] 



Statistical uncertainties in the values of the BL frac- 
tion are small and can be neglected, as confirmed by the 
dispersion of the points in FigEl More important are the 
systematic ones associated with the imprecise knowledge 
of the shape of the BL spectrum and its possible varia- 
tions along the Z-track. In order to probe the amplitude of 
these uncertainties on the horizontal branch we explored 
the range of la errors of the parameters of the comptt 
model obtained in the section 3. The associated change of 
the boundary layer fraction did not exceed w 0.05 — 0.07 
in the units of Fig. As for the normal branch, there 
is some difference between the frequency resolved spec- 
tra of GX340-I-0 on the horizontal and normal branch, 
the latter being better described by the Wien spectrum 
with kT=2.4 keV than by the spectrum of saturated 
Comptonization (cf . thick and thin dashed lines in Figl^j) . 
This difference could be related to weak dependence of 
the boundary layer spectrum on the global mass accretion 
rate (Gilfanov. Revnivtsev & Molkov 200^. Substitution 
of the comptt model by the blackbody or Wien spectrum 
with kT = 2.4 keV, decreases the boundary layer fraction 
by « 0.1. As it is clear from these numbers, systematic 
uncertainties might change the FigEl in details but does 
not affect the general trend in a significant way. 

The Fig01 and suggest that the boundary layer frac- 
tion decreases along the Z-track and it is smaller on the 
normal branch that on the horizontal branch. As demon- 
strated above, this conclusion is rather robust, as long 
as the assumption regarding the constancy of the bound- 
ary layer spectrum is approximately correct. Even if the 
diskbb model is not applicable on the normal branch any 
more, the ^ 5-fold decrease of the £^ ^ 15 keV flux be- 
tween Sz = and 5*^ = 2 effectively constrains the bound- 
ary layer contribution on the normal branch (FigEj). As 
the variability at / ^ 1 Hz is primarily associated with 
the boundary layer emission, the decrease of the boundary 
layer fraction along the Z-track also explains well-known 
decrease of the level of aperiodic and quasi-periodic vari- 
ability. 

Although no complete physical interpretation of the 
observed behavior can be offered, we mention several pos- 
sibilities. One of these is that the general structure of the 
accretion flow does not change signiflcantly and ~ 50% of 
the energy is always released on, or very close to the neu- 
tron star surface. The apparent decrease of the boundary 
layer fraction on the normal branch is a result of its ge- 
ometrical obscuration by, for example, the geometrically 
thickened accretion disk. An alternative possibility is that 
at high values of the mass accretion rate M ~ MEdd a sig- 
nificant modification of the accretion flow structure might 
occur and its division into two geometrically distinct parts 
- boundary layer and accretion disk, might become inap- 
plicable. Namely, due to non-negligible pressure effects the 
deceleration of the orbital motion of the accreting matter 
from Keplerian frequency to the neutron star spin fre- 
quency would take place in a geometrically extended re- 
gion with the radial extend of AR ~ i?NS- In this case, the 
observed decrease of the boundary layer fraction could re- 
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fleet actual decrease of the fraction of the energy released 
on the neutron star surface with the rest of the energy be- 
ing released in the extended transition region. Finally, it is 
also possible that the apparent decrease of the boundary 
layer fraction as well as the peculiar behavior of the ac- 
cretion disk temperature (FigEJ is the consequence of the 
complete break down of the model on the normal branch 
of the color-color diagram. 

Sco X-1 and GX17-f 2 appear to have higher boundary 
layer fraction at the end of the normal branch than other 
three sources. Another indication of the possibly higher 
contribution of the BL emission at 5*2 ~ 2 in Sco X-1 
and GX17-I-2 is that they typically show stronger fast 
variability at the normal and flaring branch than Cyg X-2, 
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be mentioned though, that the appearance of the Fig^lis 
a direct consequence of the definition of the Sz coordi- 
nate, in particular the definition of the ~ 2 reference 
point as the turning point from the normal to the flaring 
branch. If, for example, the BL fraction would be plotted 
against the hard color, all five sources will be almost 
indistinguishable . 

5.3. Shape of Z-diagram 

Motivated by these results we use the two-component 
spectral model introduced in the beginning of this section 
to explain behavior of the Z-sources on the color-color dia- 
gram. Unlike before, for the accretion disk spectral compo- 
nent we adopt here the general relativistic accretion disk 
model by Ebisawa et al. (1991) (the grad model in the 
XSPEC fitting package) which explicitly includes depen- 
dence on the mass accretion rate. The mass of the com- 
pact object and binary system inclination were fixed at 
-^NS = 1.4M0 and i = 60°. Using this two-component 
spectral model we can calculate the position in the color- 
color diagram as a function of the mass accretion rate and 
the boundary layer fraction. This simple model is an at- 
tempt to understand general tendencies of the Z-diagram, 
rather than to construct its precise quantitative descrip- 
tion. In particular, the fact is ignored, that the relation 
between the accretion disk temperature and its flux is 
more complex than predicted by the standard model of 
the geometrically thin disk (FiglSJ. 

The results are presented in Figd The two dashed 
curves in the figure show the evolution of spectral col- 
ors with the increase of the mass accretion rate of the 
disk component for two different values of the boundary 
layer fraction. The upper curve assumes the BL-to-disk 



flux ratio of 0.8, in the lower curve the BL fraction equals 
zero. On each curve positions where the mass accretion 
rate equals M = 0.5, 1, 2, 3, 4 x 10^^ g/s are marked. The 
evolution of colors assuming a rapid change of the BL con- 
tribution from i^BL/-Fdisk = 0.8 to zero approximately at 
M = 2 X 10^^ g/s is shown by thick solid line. 

The general shape of the Z-track can be reproduced in 
the model as a result of variation of two parameters - the 
mass accretion rate and BL fraction. The mass accretion 
rate increases along the Z-track. The Z-shape of the track 
is defined by the variation of the BL fraction, which de- 
creases along the normal branch from the value of ~ 50% 
expected in the "standard" theories to a small number of 
the order of '^zero at the end of the normal branch. 

Quantitatively, the exact value of the Al, correspond- 
ing to the transition from the horizontal to the normal 
branch depends on the disk model (grad) parameters - 
the binary system inclination, the mass of the neutron 
star and the spectral hardening factor. For our choice of 
parameters, it equals M ~ 2 • 10^^ g/sec, i.e. is of the or- 
der of the Eddington critical value for a 1.4Mq neutron 
star, when we can expect that the structure of the in- 
ner accretion flow can start to change. It is interesting to 
note that similar conclusions (decrease of the observable 
NS emission due to possible inflation of the inner disk) 
were obtained from studies of simultaneous observations 
of Z-source s in radio, optical UV and X-ray ranges (e.g. 
iHasinger et al. 1990; Vrtilek et a l. 1990) 

The location of the horizontal branch on the color- 
color diagram depends on the boundary layer fraction at 
smaller M. Theoretically, it could be used to determine 
the latter from observations which would be of great help 
for the theory. In particular, our spectral model suggests 
Fbl/Fad ~ 0.8. In practice, however, the direct compar- 
ison of this number with theoretical predictions for the 
energy release in the disk and boundary layer is compli- 
cated, because the observed flux ratio Fbl / -Fdisk is modi- 
fied by geometrical factors and anisotropy of the AD and 
BL emission diagrams. 



6. Summary 

Using the archival data of RXTE observations of several 
bright LMXBs we study the boundary layer emission and 
the structure of the accretion flow in this sources. Our 
sample includes 2 atoll sources in the soft spectral state 
and 5 (all but one) Z-sources. Their luminosities range 
from ^ 10'^^ erg/s (atoll sources) to ^ 2 • 10'^* erg/s (Z- 
sources). Our results are summarized below. 

— We construct Fourier-frequency resolved spectra (en- 
ergy spectra of variable emission) at / ^ 1 — 5 Hz 
for five sources, which data combines sufficient timing 
and spectral resolution. Within the statistical uncer- 
tainties all the spectra have the same shape (FigEJ. 
This is a remarkable result as the average spectra of 
the atoll and Z-sources differ significantly and their lu- 
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Fig. 7. The color-color diagram of 5 Z-sources. Overlayed 
on the data are the tracks predicted by the model (section 
I5.3|l . The thick two dashed lines show evolution of the 
colors with change of the accretion rate in the disk for 
two different values of the BL fraction: 44% (upper) and 
zero (lower) The values of the mass accretion rates M of 
the disk component are marked in units of 10^* g/s. The 
thick solid line shows the Z track with transition mass 
accretion rate 2 x 10^^ g/s. 



minosities range from ~ 10'^^ erg/s (atoll sources) to 
~ 2 ■ 10^* erg/s (Z-sources). 

This result fits in the pi ct ure suggested by 
iGilfanov. Revnivtsev fc Molkovl l)2003^ and briefly 
outlined in the Introduction. This picture implies, 
in particular, that the frequency resolved spectrum 
equals the energy spectrum of the boundary layer emis- 
sion. The following is based on this assumption. 

— The boundary layer spectrum in atoll sources and on 
the horizontal branch of Z-sources (i.e. at presumably 
sub-Eddington accretion rates) can be approximated 
by the spectrum of saturated Comptonization (section 
O. Its high energy {E >10 keV) cut-off corresponds to 
the black body spectrum with temperature of kT ^ 2.4 
keV (FigEJ. 

— Assuming that the boundary layer is radiation pressure 
supported we constrain the neutron star mass and ra- 
dius. For a 1.4M0 neutron star and the spectral hard- 
ening factor of 1.6-1.8 the NS radii are in the range of 
i?NS 9 - 14 km (FigEJ. 

— We attempt to relate the motion of Z-sources along 
the Z-track to changes in the values of the physically 
meaningful parameters. Our results tentatively suggest 



that the contribution of the boundary layer component 
to the observed emission decreases along the Z-track 
from the conventional value of ~ 50% on the horizontal 
branch to a rather small number at the end of the nor- 
mal branch. The main difference of our approach from 
previous attempts is in the apriori knowledge of the 
shape of the boundary layer spectrum. This allowed 
us to avoid ambiguity of the spectral decomposition 
into boundary layer and disk components which was 
one of the major problem in previous studies. 
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